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ABSTRACT 

We report the first detection of asymmetry in a supernova (SN) photosphere based on SN light 
echo (LE) spectra of Cas A from the different perspectives of dust concentrations on its LE ellipsoid. 
New LEs are reported based on difference images, and optical spectra of these LEs are analyzed and 
compared. After properly accounting for the effects of finite dust-filament extent and inclination, we 
find one field where the He I A5876 and Ha features are blueshifted by an additional ^4000 kms~^ 
relative to other spectra and to the spectra of the Type lib SN 1993J. That same direction does 
not show any shift relative to other Cas A LE spectra in the Ca II near-infrared triplet feature. We 
compare the perspectives of the Cas A LE dust concentrations with recent three-dimensional modeling 
of the SN remnant (SNR) and note that the location having the blueshifted He I and Ha features 
is roughly in the direction of an Fe-rich outflow and in the opposite direction of the motion of the 
compact object at the center of the SNR. We conclude that Cas A was an intrinsically asymmetric SN. 
Future LE spectroscopy of this object, and of other historical SNe, will provide additional insight into 
the connection of explosion mechanism to SN to SNR, as well as give crucial observational evidence 
regarding how stars explode. 

Subject headings: ISM: individual (Cas A) — supernovae: general — supernova remnants 
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1. INTRODUCTION 

Light echoes (LEs) are the scattered light of a transient 
event that reflects off dust in the interstellar medium. 
The extra path length of the two-segment trajectory re- 
sults in LE light arriving at an observer significantly later 
than the undelayed photons. Such circumstances provide 
exciting scientific opportunities that are extremely rare 
in astronomy — specifically, to observe historical events 
with modern instrumentation and to examine the same 
event from different lines of sight (LoS). We have previ- 
ously employed LEs to take advantage of the time delay, 
identifying systems of such LEs associated with several- 
hundred-year old supernova (SN) remnants (SNRs) in 
the Large Magellanic Cloud and subsequen tly spectro- 
scopically cla ssifying the super novae (SNe; Rest et al. 
2005b, 2008a). Our work represented the first time that 
the spectral classification of SN light was definitively 
linked with a S NR. The analysis of th e X-ray spectrum 
of that SNR by [Badenes et all ( 12008D provided confir- 
mation of both its classification (Type la) and subclass 
(high luminosity). Clearly, LE spectroscopy is a power- 
ful technique for understanding the nature of SNe in the 
Milky Way and Local Group galaxies. 

Cas A is the b rightest extrasolar radio source in the sky 
(Ryle & Smith 1948| ) and the younge st (age ^330 yr) 
Milky Way core-collapse SNR (Stephenson & Green 
2002). Its distance is approximately 3.4 kpc. Dynami- 



probe the SN from significantly different directions. The 
onl y time this te c hniqu e has b een applied previously was 



cal measurements of the S NR indicate that th e explosion 
occurred in year 1681 it 19 dFesen et al.|2006aD ; we adopt 
this date for age calculations in this paper. A single 
historical account of a sighting in 1680 by Fl amsteed is 
attributed to the Cas A SN (Ashworth 1980), although 
this has been disputed (Kamper 198 0). Cas A is the 
youngest of the certain historical CC SNe and is thus an 
e xcellent target for L E studies. 



Krause et al. (2005) identified a few moving Cas A fea- 
tiires (called ''infrared echoes") using infrared (IR) im- 
ages from the Spitzer Space Telescope, the result of dust 
absorbing the SN light, warming up, and reradiating light 
at longer wavelengths. Their main scientific conclusion, 
that most if not all of these IR echoes were caused by 
a series of recent X-ray outbursts from the compact ob- 
ject in the Cas A SNR, was incorrect because they did 
not take into account that the apparent motion strongly 
depen ds on the inclination of the scattering dust fila- 
ment dPwek fc Arendt||2008| [Rest et al.||2011| ). Rather, 
the echoes must have been generated by an intense and 
short burst of ultraviolet (UV) radiation associated with 
the breakout of the SN s hock through the surface of the 
progenitor star of Cas A (iDw ek fc Arendt 2008) . 

The first scattered LEs of Galactic !SiNe associated with 
Tycho ' s SN and the Cas A SN were disco vered by [Rest 



?t al. l I 

mm 



et al.| ( 2007 2008b ) . Contemporaneously, [Krause et al. 



obtained a spectrum of a scattered optical light 
echo spatially coincident with one of the Cas A IR echoes, 
and identified the Cas A SN to be of Type lib from its 
similarity to the spectrum of SN 1993 J, the very well- 
obser ved and prototypical example of the SN lib class 
(e.g.,|Filip penko et al."1993^ 'Richmond et al.||1994[ [FiT 
ippenkoTe t al. 1994; Matheson et al. 2000). 

'I'he discovery and spectroscopic follow-up observations 
of different LEs from the same SN allow us to benefit 
from the unique advantages of LEs — their ability to 



by [Smith et al.| ([2001] [2003| , who used spectra of the 
reflection nebula of 77 Cannae to observe its central star 
from different directions. Dust concentrations scattering 
SN light lie at numerous different position angles and at 
different radial distances from the observer. This is il- 
lustrated in Figure [1] which shows the Cas A SNR (red 
dots), the scattering dust (brown dots), and the light 
paths. We denote the light echoes as LE2116, LE2521, 
and LE3923, where the number represents the ID of the 
grid tile of the search area in this region of the sky. 
LE2521 and LE39 23 are newly discov ered and LE2116 
was discovered by Rest et al. (2008b). The scattering 
dust of LE3923 is more than 2000 ly m front of Cas A, 
much farther than any other scattering dust; thus, we 
show only part of its light path. 

The study of a single SN from different LoS is partic- 
ularly relevant for the Cas A SNR. With observations 
of the LEs shown in Figure [l] we are positioned to di- 
rectly measure the symmetry of a core-collapse SN, and 
compare it to the structure of the remnant. The ther- 
mal X-ray emission as well as the optical emission of the 
Cas A SNR is very inhomogeneous, with large Fe-rich 
and Si-rich outflows being spatially disti nct, indicating 



2000[ [Hwang e t al.|2004l|Fesen et al.|2006a| rT: 

of these outflows and their individual relevance 



that the SN explosion was asymmetric ([Hughes et al. 

he nature 
to the SN 

explosion is still debate d (Burrows et al. 2005; Wheeler] 



et al 



2008 DeLaney et al. 2010). Tananbaum (1999 



)99J 



detected in Chandra images a compact X-ray source' 
from the Cas A SNR center; it is an excellent candidate 
for b eing the neutron star produced by the SN explo- 
sion (Tan anbaum[[1999| [Fesen et al][2006b[ ). The po- 
sition angle of the X-ray source is oft by only ^ 30° 
from the position angle o f the southeast Fe-rich struc- 
ture (Wheeler et al.|2008 ), and has a projected apparent 
motion of 350 kms . In ^of this paper, we report two 
new LE complexes associated with Cas A, and we show 
optical spectra of these LEs as well as one of the previ- 
ously known LEs. With these data we are able to view 
three distinct directions, where each dust concentration 
probes different hemispheres of the SN photosphere. In 
Rest et al.[ ([2010") we have introduced an innovative tech- 



nique for modeling both astrophysical (dust inclination, 
scattering, and reddening) and observational (seeing and 
slit width) effects to measure the light-curve weighted 
window function which is the determining factor for the 
relative time-weighting of the observed (integrated) LE 
spectrum. We apply this technique in ^to similar well- 
observed SNe to produce appropriate comparison spec- 
tra. Such spectra are necessary to compare the Cas A 
spectra to other SNe as well as to compare the LE spec- 
tra to each other. We show that the Cas A SN was in- 
deed very simil ar to the proto t ypical Type lib SN 1993 J, 



as claimed by Krause et al. 
strate that, despite the excel 



(2008|). In Pwe demon- 
ent agreementbetween the 



Cas A spectra and SN 1993J, one LE has a systematically 
higher ejecta velocity than either SN 1993J or the other 
LEs, revealing that Cas A was an intrinsically asymmet- 
ric explosion; observers from different directions would 
have viewed a "different" SN spectroscopically. In this 
section we also discuss the implications of our finding for 
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Fig. 1. — North is toward the positive-y axis (up), east is toward the negative-x axis (left), and the positive- 2; axis points toward the 
observer with the origin at the SNR. The SN to LE cone apex distance, in hght years, is half the interval of time since the SN explosion. 
The red and brown circles indicate the SN and scattering dust, respectively. The black lines show the path of the light scattering from the 
LE-producing dust concentrations. The scattering dust of LE3923 is more than 2000 ly in front of Cas A, much farther than any other 
scattering dust; thus, we show only part of the light path, and do not include its dust location. The top-left panel shows an Chandra X-ray 
image (Hwang et al. 2004), with the projected light path from SN to scattering dust overplotted (gray arrows). The red arrow indicates the 
X-ray compact object and its apparent motion. In this false-color image, red corresponds to low-energy X-rays around the Fe L complex 
(~1 keV and below), green to mid-energy X-rays around the Si K blend (~2 keV), and blue to high-energy X-rays in the 4-6 keV continuum 
band between the Ca K and Fe K blends. 
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both other historical SNe and for core-cohapse SNe and 
their explosions. 

2. OBSERVATIONS Sz REDUCTIONS 
2.1. Imaging 

We reobserved LEs discovered in a campaign of several 
observing runs on the Mayall 4 m telescope at Kitt Peak 
National Ob servatory (KPNO) starting in 2006 (Rest 



et al.||2008 bl). As described in that paper, the Mosaic 
imager, which operates at the //3.1 prime focus at an 
effective focal ratio of //2.9, was used with the Bernstein 
VR broad-band filter (kl040) which has a central wave- 
length of 5945 A and a full width at half-maximum in- 
tensity (FWHM) of 2120 A. The images were kernel- and 
flux-matched, aligned , subtracted , and masked using the 



naitis et al. 2007)7 



SMSN pip ehne (Rest et al.|2005a 

m 



Garg et 



LE211D had been previously 



al.|2007l 

dously di 



Mik- 



ered and reported (Rest et al. 2008b), whereas LE2521 
and LE3923 are light echoes discovered on 2009 Septem- 
ber 14 and 16, respectively (UT dates are used through- 
out this paper). 

Figure [2] shows the LEs and the adopted spectroscopic 
slit positions. For the LE2116 mask design, we used 
images from KPNO obtained on 2009 September 14, 
~1 week before the spectroscopy was obtained. Note 
that we had a third slit in the high surface brightness re- 
gion between slits A and B. Unfortunately, the LE filled 
the entire slit, making proper sky subtraction impossi- 
ble. Sky subtraction using the sky observed in other slits 
produced inaccurate spectra with large residuals near the 
sky lines. As a result, we have not included the spectra 
from this slit in our analysis. For LE3923, we used im- 
ages obtained at the Apache Point Observatory 3.5 m 
telescope on 2009 October 16 with the SPIcanQ CCD 
imager in the SDSS-r' filter and processed in a similar 
manner with the SMSN pipeline. Since the apparent 
motion of the LEs is ~30'' per year, we measured the ap- 
parent motion and adjusted the slit position accordingly. 
These position adjustments amounted to . LE2521 
was sufficiently compact and bright that long-slit spec- 
troscopy was favored over masks. A summary of the 
spectroscopic and geometric parameters is given in Ta- 
bles [1] and |2j respectively. For the light-echo profile fits 
described in ^ we used the deep images in good seeing 
from 2009 September 14-16. 

2.2. Spectroscopy 

We obtained spectra of the LEs with the Low Re so- 
lution Imaging Spectrometer (LRIS; Oke et aL]|1995 ) on 
the 10 m Kec k I telescope and wi th the Blue Channel 
spectrograph QS chmidt et al.||1989| ) on the 6.5 m MMT 
(see TablelT]). For the LRIS observations, slit masks were 
designed to maximize the efficiency of the telescope time. 
Standard CCD proces sing and spectrum extraction were 
performed with IRAFp^ The data were extracted us- 



ing the optimal algorithm of Horne (1986). Low-order 
polynomial fits to calibration-lamp spectra were used to 

http: / /www. apo.nmsu.edu /arc35m/Instruments/SPIC AM / . 

IRAF: the Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astron- 
omy (AURA), Inc., under cooperative agreement with the National 
Science Foundation (NSF). 
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Fig. 2. — Image stamps of the LE discovery images and slit posi- 
tions for spectroscopy for the three LE complexes LE2116, LE2521, 
and LE3923 (along horizontal panels). Left vertical panels are the 
reference template image at an earlier epoch (on the date shown 
in the box at top). Middle vertical panels are the image fields at 
a later epoch, centered here on the LE position. Right vertical 
panels are the difference images of the respective recent image and 
the earlier reference template to isolate the transient fiux of the 
LE. For all image stamps north is up and east is to the left. The 
position and size of the slit used for spectroscopy is indicated by 
the rectangular overlay. 



establish the wavelength scale, and small adjustments 
derived from night-sky lines in the object frames were 
applied. We employed our own IDL routines to flux cal- 
ibrate the data and remove telluric lines using the well- 
expose d continu a of s pectrophotometric standard stars 



(Wade fc Horn e 1988; Matheson et aI]|2QQQ[ [F oley et al. 
2003). For LE2116, we combine slits AT^and B usmg 



weights of 0.9 and 0.1, respectively. We present our LE 
spectra in Q 

3. GENERATING COMPARISON SPECTRA 

We have recently determined that modeling the dust- 
filament properties (e.g., the dust width) is re quired to 
accurately model LE spectra (Rest et al. 2010). All pre- 
vious studies have neglected this aspect. Because the 
dust filament that reflects the LE has a nonzero size, it 
translates into a window function over the time domain 
as the SN light traverses over it. If this window func- 
tion is narrower than the SN light curve, the observed 
LE spectrum is affected. We find that the observed LE 
spectrum is a function of dust-filament thickness, dust- 
filament inclination, seeing, spectrograph slit width, and 
slit rotation/position with respect to the LE. All of the 
noninstrumental parameters can be determined with an 
analysis of the LE in images at different epochs, and an 
individually modeled spectral template for a given ob- 
served LE can be constructed a priori wi t hout using the 
actual observed LE spectrum. Rest et al. ( 2010[ ) describe 
this process in detail, showing that the observed LE spec- 
trum is the integration of the individual spectra weighted 
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TABLE 1 







R.A. 


deck 


PA SNR-LE 




Seeing^ 




Width^ 


Length^ 






LE 


Telescope 


(J2000) 


(J2000) 


n 


UT Date^ 


n 


ID^ 






n 


UT Mask^ 


2116 


Keck 


23:02:27.10 


+56:54:23.4 


237.79 


20090922 


0.81 


A 


1.5 


2.69 


0.0 


20090914 


2116 


Keck 


23:02:27.67 


+56:54:07.7 


237.79 


20090922 


0.81 


B 


1.5 


3.59 


0.0 


20090914 


2521 


Keck 


23:12:03.86 


+59:34:59.3 


299.10 


20090922 


0.79 




1.5 


4.32 


41.0 




2521 


MMT 


23:12:03.88 


+59:34:59.2 


299.10 


20090921 


0.70 




1.0 


4.20 


30.0 




3923 


Keck 


00:21:04.97 


+61:15:37.3 


65.08 


20091023 


0.89 




1.5 


7.72 


0.0 


20091016 



^ Parameters for spectroscopic slit. 

^ UT date of image used for mask design. 



TABLE 2 





5^ 










LE 


n 


(ly) 


(ly) 


(ly) 


(ly) 


2116 


3.376 


628 


531 


-335 


437 


2521 


1.642 


318 


278 


155 


-10 


3923 


7.587 


1204 


-1092 


507 


2045 



^ Angular distance between SNR and LE. 

^ Distance between SNR and LE in the plane of 

the sky (p^ = x'^ -\- y'^). 

*^ Coordinates x,y, and z of the LE with origin 
at the SNR. The positive x and y axis are in the 
plane of the sky toward West and North, respec- 
tively. The positive z axis is along the line of 
sight from the SNR toward the observer. We as- 
sume a distance to the Cas A SNR of 3.4 kpc and 
1681 AD as the time of explosion. 

with an effective light curve, which is the product of 
the hght curve with the dust-filament-dependent window 
function. We emphasize that these window functions are 
different for every LE location. The top-left panel of Fig- 
ure |3] shows the window functions for the LEs for which 
we have spectra using the SN 1993 J spectral library ( Jef-j 



J: ^, - . 

fery et al.|1994[ |Barbon et al.il 995; Fr ansson et al.|2i 

and light 'curve (|KichmQnd et al.^ 1995] ^ We refer the 



reader to Rest et"al.| ( |2U1U[ ) for a detailed description of 
how the ettective lignt curves are derived. Examples of 
the fitted parameters are given in Tables 1 and 2 and 
shown in Figure 12 of Rest et al. (2010) 

Even though LE2521 and LE2116 have comparable 
dust widths, LE2521 has a considerably wider window 
function than LE2116 since its dust-filament inclination 
[a = 54°) is significantly larger than that of LE2116 
{a = 9°), causing the projected SN light-curve shape to 
be "squashed." LE3923 has an inclination similar to that 
of LE2116, but its scattering dust is thicker, leading to a 
wider window function. We note that the dust-filament 
inclination and width are the two most influential pa- 
rameters affecting the window function. 

The middle-left and bottom-left panels of Figure 
show the effective light curves of SN 1993 J an( 
SN 2003bg, respectively. We constructed the effective 
lig ht curve of SN 2 Q03bg using the spectra and light curve 
in Hamuy et al. (2009|). In general, the rise time of a 
SNf 



iamuy et al.| ( |2QQ9[) . 

is shorter than 50 days. Since the window functions 
have, in most cases, FWHM > 30 days (see the top- 
left panel in Figure |3|, the rising wing of the effective 
light curves and the real light curve are similar (see the 
middle-left and bottom- left panels in Figure [3|, and con- 
sequently the impact on the integrated spectra is small. 
However, the decline time of a SN is significantly longer, 
and therefore the window function has a much more pro- 



found impact at late phases: it cuts off the contribu- 
tion of spectra with phases > 100 days. Even though 
the intrinsic brightness of the SN at late phases is much 
fainter than at peak, these late-phase spectra still sig- 
nificantly contribute to the integrated spectra since they 
are completely dominated by a few persistent lines. This 
is illustrated in the top-right and bottom-right panels of 
Figure [3j which show the integrated spectra of SN 1993 J 
and SN 2003bg, respectively. The late-time spectra of 
SNe lib such as SN 1993J and SN 2003bg are dominated 
by [O I] AA6300, 6363, [Ca II] AA7291, 7324, and the 
Ca II NIR triplet. Note that the black spectra, which 
are weighted by the original, unmodified light curve, are 
much stronger at these lines compared to the other spec- 
tra, which are integrated using the effective light curves. 
This figure clearly illustrates the impact and importance 
of the effective light curve on the contribution of the late- 
phase spectra to the observed integrated spectrum. 

3.1. Integrating and Fitting the Spectral Templates to 
the Observed LE Spectra 

We numerically integrate the spectral t emplate s of 
SN 1993J (|Jeffery et al. 199l[ | Barbon et al. 1^1995 j RidT] 



mond et aTTT 



Hamuy et a. 



mo 

hi 



Fransson et al.||20 05|), ^5]T^(MEg 



2009[ ), and SN 2008ar lUhornock et al. 
010a| ) by weignting t he individual epoch s with the light 



curve as described bylRest et al. (2008a!), with the only 
difference being that the real lignt curve is replaced with 
an effective lig ht cu r ve as described in ^ and in more 



detail by Rest et al. (2010) 



Several authors have addre ssed the single - scattering 



approxima tion for LEs 
1986 



Couderc'^T939! "Chevalier 
DS9; Sugerma^nmr 



Emmering & Chevalier 
Patat 2005 )'^ The total surface brightness depends on 
parameters such as the intrinsic brightness of the event, 
dust density and thickness, and others. For our purposes, 
we only need to consider the parameters that attenu- 
ate the spectrum, which are (1) the forward scattering 
described by the integrated scattering function S{X,0) 
as derived by Sugerman (2003), and (2) the reddening 
due to the extinction by Galactic dust. The attenua- 
tion by forward scattering can be determined a priori 
using the observed angular separation of the LE from 
Cas A. The only uncertainties introduced here are due 
to the assumed age and distance for Cas A. The redden- 
ing cannot be determined independently and therefore 
we fit for it. We fit the template spectra to the observed 
LE spectra with the only free parameters being the nor- 
malization and the reddening. Figure [4] compares the LE 
spectrum of LE2521 (blue line) to various SNe lb, Ic, and 
II using the window function for LE2521. SN 1993 J and 
SN 2003bg quite clearly provide the best fit. 
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Fig. 3. — Red, cyan, and blue indicate LEs LE2521, LE3923, and LE2116, respectively. Top left: window function for each LE. Middle 
left: effective light curves of SN 1993J for each LE. The unmodified light curve of SN 1993J is shown in black. Epochs with spectra 
are indicated by the dashed lines. Bottom left: same as middle left, but for SN 2003bg. Top right: integrated SN 1993J spectra, where 
the integration is weighted by the respective effective light curve. Note that the black spectrum is weighted by the original, unmodified 
SN 1993J light curve. The differences in the spectra are mainly due to the strong lines of [O l] AA6300, 6363, [Ca II] AA7291, 7324, and the 
Ca II NIR triplet in the late-phase spectra. Bottom right: same as top right, but for SN 2003bg. 



4. COMPARISON OF THE CAS A LIGHT ECHOES 

In Figure |5) we present four spectra of Cas A LEs, 
corresponding to three different directions. LE2521 was 
much brighter than the other LEs, and as a result, its 
spectrum has a higher signal-to-noise ratio (S/N) than 
the others. The spectrum of LE3923 has the lowest S/N, 
but it still displays the same spectral features as the other 
spectra. All three directions have very similar spectra 
with the main differences being the continuum shape, 
which we have shown to be a result of different scattering 
and reddening for the various directions. The line profiles 
do exhibit some differences, but as shown in ^ different 
window functions can create significant differences in the 
spectra. 

To properly compare the LEs from different directions, 
one must make a comparison to another SN, which re- 
moves the effect of the window function. Assuming that 
the window functions are correct, if a LE spectrum is the 
same as the integrated spectrum of a SN, then the two 
objects likely looked the same at maximum light (there is 
obviously some degeneracy between the different epochs 
within the window function). If the two spectra are dif- 
ferent, then either their maximum-light spectra or light 
curves were likely different; regardless, it is a direct in- 
dication that the two objects appeared different along 
those LoS. Moreover, if one LE matches the integrated 
spectrum of a particular spectrum, but a LE from a dif- 
ferent direction does not, then observers along the two 
LoS would have seen "different" SNe. 

For Cas A, we have compared the LEs to spectra of 
SNe 1993 J, 2003bg, and 2008ax (see FigM. These are 
the best-observed SNe lib for which we had light curves 
and a good time series of spectra. As described above, 
SN 1993 J is an iconic SN lib that was extensively ob- 



served (e.^., Filippenko et al.*1993URichmond et al.|1994 
|l^ ilippenko et al. 1994; Matheson et al. 2000) and whose 
progenitor system was identified in pre-explosion images 
(Podsiadlowski et al. 1993 ; Aldering et al. 1994) , and its 
bmary co mpanion was ident ified through late-time spec- 
troscopy (i Maund et al.|2004 ). SN 2003bg was a luminous 
event (both op tically and in the rad io) with broad lines 



at early times (Soderberg et al. 12006 Hamuy et al.|200^ 
Mazzali et al. 2009). SN 2008ax was discovered hour^ 



)). SN 2008ax was discovered hours 
alter explosion , which enabled significant follow-up ob- 



servations 



Pastorello et al. |2008[ [Roming et al.||2009 
2010b ) 71ncludin^ extensive spectropo- 
Chornock et al. 2010b) indicating that the 

Fur- 



Chornock et al 
larimetry 

outer layers of the eject a were rather aspherical. 
thermore, its progenitor system was also detected in pre- 
explosion images (Crockett et al. 2008). 

Examining Figure [5[ it is clear that the LE spectra 
in all directions are similar to each other. In addition, 
the LE spectra are similar to spectra of SNe 1993 J and 
2003bg, indicating that no dramatic differences are seen 
from the various directions (such as a different spectral 
classification). SN 2008ax is a worse match to the LEs, 
having significantly weaker Ha emission than the other 
objects. Since our spectral seri es of SN 2008ax (both 
CfA and published spectra from Chornock et al. 2010b) 
is not as exhaustive as that of SNe 1993J and z003bg, 
we do not know if the differences are intrinsic or the 
result of incomplete data. Because of the poor match 
and the differences not necessarily being physical, as well 
as for clarity, we do not show in Figure [5] the comparison 
spectra generated from SN 2008ax. 

Despite all LEs being quite similar, a detailed com- 
parison of the LE line profiles to those of SN 1993J 
in Figure [g] indicates that while two directions (LE2116 
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Fig. 4. — The LE spectrum of LE2521 (blue line) compared to 
various SNe lb, Ic, and II using the window function shown in 
Figure [3] SN 1993 J and SN 2003bg unambiguously give the best 
fit. 



and LE3923) are consistent with the hne velocities of 
SN 1993J, one direction (LE2521) has higher-velocity 
absorption for the He I A5876 and Ha features. Specif- 
ically, the minima of the features are blueshifted by an 
additional 4000 and 3000 kms~^, respectively. In ad- 
dition to the blueshifted absorption minima, the entire 
profile of the features seems to be shifted blueward, in- 
cluding the emission component of the P-Cygni profile. 
The Ca II NIR triplet may have a slightly larger velocity 
for LE2521 than in the other directions, but this does 
not appear to be significant. We estimate that the un- 
certainties in these lines are small and not significantly 
more than ~100 kms~^. 

Interestingly, the comparison spectra from SN 2003bg 
show Ha velocities that are larger than those of Cas A 
for LE2116 and LE3923, but the velocity matches that of 
Cas A for LE2521. For this direction, the profile shape 
is also matched. The LE2521 He I velocity (and full 
profile) appears to be larger than that of SN 2003bg, 
which has a similar velocity to that of SN 1993 J. This 
indicates that although SN 2003bg is an excellent match 
to the Ha feature of LE2521, there are still differences, 
and SN 2003bg is not a perfect analog. SN 2003bg has 
a slightly higher velocity for the Ca II NIR triplet than 
SN 1993 J, which makes it more consistent with LE2521, 
but we reiterate that the difference is not significant. 



5. COMPARISON OF THE LIGHT ECHOES TO THE SNR 
STRUCTURE 

The structure of the Cas A SNR is very complex — 
from the large scales, which show spatially distinct Fe- 
rich and Si-rich regions, to the small scales, which are 
dominated by many filaments and bright knot-like fea- 
tures in the optical and X-ray bands. Overall, the ejecta 
emission from the SNR appears very turbulent and in- 
homogeneous, similar to other remnants of core-collapse 
SNe (e .g., G292.0+1.8, Park et al.| |2007| ). [Hughes et al] 
(2000) argue that the ejecta have overturned during the 
explosion, since the Fe-rich material is ahead o f the Si- 
rich materia l in the southeast (SE; however, see DeLaney 



et al. 



2010 ). The Si-rich regions in the northeast (NE) 

an d southwest (SW) sh ow structure that was interpreted 
by Hwang et al | (|2004|) as indicative of a bipolar j et sys- 



teni. I'hese ''jet-like'' outflows were also seen by Fesen| 
et al.| ( [2Q06a ) in the Hubble Space Telescope (HSTJ cen- 
sus or last-moving optical knots of the intermediate- mass 
elements N II and O II, further supporting this NE-SW 
bipolar jet interpretation. 

While this "jet-counterjet" structure is clearly a direc- 
tional outflow associated with the ejecta (and not the 
result of some pre-ex isting bipolar cavity in the circum- 



stellar medium; e.g., Blondin et al. 



1996| ), it seems that 



the outflow was not energetic en ough to power the entire 
explosion (Laming et al. 2006). In addition, the pro 



jecte d motion o f the co mpact object detect ed in the X- 
rays ( Tananbaum|[l999 Fesen et al. 2006b| ), thought to 
be the neutron star or black hole created by the explo- 



y tne expL 

sion, is perpendicular to the NE-SW outflows. Most jet- 
powered explosion models predict that the neutron-star 
kick would be roughly aligned with the jet axis. Fur- 
thermore, the Fe-rich outflows are concen trated int o two 
regions in the SE and NW ([Hughes et al. ^ ,200Q] iWill- 



ingale et al. 



flow blueshi: 



2002[ [Dewey et al.||2007l ), with the ^E_Qut^ 
^tea~and the INW outflow redshifted fWill- 



ingale et al. 2002). This has given rise to a new picture 
that there is an axis in the Cas A SNR from the SE to 
the NW at a position angle of ~125° aligned with the 
Fe-rich knots. In this model, the Si-rich structures to 
the NE and SW are only secondary features caus ed by 
instability-power ed flows from an equ atorial torus ([Bur- 
rows et al.J|2005l iWheeler et al.""2008 



0Q8|. 

In a recent study, DeLaney et al. (5oTo| used the tech- 
nique of Doppler imaging to give the most complete and 
updated three-dimensional (3D) model. This model is 
derived from IR (Spitzer), optical (ground-based and 
HST)^ and X-ray {Chandra) data. |DeLaney et al. (2010) 
used velocity information from specific spectral features 
to deproject the structure of the SNR along the radial di- 
rection perpendicular to the plane of the sky (see their §2, 
§3, and §4 for details of the data reduction and deprojec- 
tion techniques). Their work reveals a very complex 3D 
structure that can be characterized by a spherical compo- 
nent, a tilted "thick disk," and multiple ejecta outflows. 
In their model, the thick disk is tilted from the plane of 
the sky at an angle of ~25° from the E-W axis and ~30° 
from the N-S axis. This thick disk contains all the ejecta 
structures, including the most prominent outflows that 
show up best in X-ray Fe K emission, but also appear 
in other datasets. Infrared emission in the [Ar II] and 
[Ne II] lines often reveals ring-like structures, which are 
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Fig. 5. — Top panel: observed LE spectra binned to 10 A pixel" ^ (black curves). The colored curves are the variance- weighted smoothed 
spectra with the red, blue, and orange curves corresponding to LEs LE2116, LE2521, and LE3923, respectively. The LE2521 spectrum 
with the shorter wavelength range is from the MMT; the remaining spectra are from Keck. Bottom panel: variance-weighted smoothed 
LE spectra (colored curves) compared to their corresponding light-curve weighted integrated SN 1993J (solid black curves) and SN 2003bg 
spectra (dashed black curves). The LE spectra are equivalent to the spectra in the top panel. The SNe 1993 J and 2003bg spectra have 
been processed to reflect the dust scattering and reddening associated with each LE (see text for details). The spectra of Cas A are very 
similar to both SNe 1993 J and 2003bg from all three directions probed; however, there are some signiflcant differences (see Q. 

somet imes seen as broken ri n^s, at the base of these out- 
flows. iDeLaney et al.|(|2010^ note that these rings appear 



at the intersection between the thick-disk structures and 
a roughly spherical reverse shock. Some of the outflows 
are bipolar, with oppositely directed flows about the ex- 
pansion center, while others are not. In particular, the 
blueshifted emission from the Fe-rich outflow in the SE 
is clearly collimated, but the corresponding redshifted 
emission in the NW is not. 



The 3D analysis by DeLaney et al.| (2QlQ ) suggests that 
the Cas A SN explosion was highly asymmetric, with 
most of the eject a flattened in a thick disk that is slightly 
tilted from the plane of the sky and no prominent struc- 
tures perpendicular to the plane of the sky. It is possible 
that this might be due to selection effects, at least to 
some extent, because limb brightening will make struc- 
tures close to the plane of the sky easier to detect, but 
given the quality of the individual datasets, it is unlikely 
that any prominent outflows along the radial direction 
would have been missed completely by [DeLaney et al.| 



([2010|. 

In Figure [7| we superimpose the LoS fro m our LEs onto 
three of the^Ar II] and Fe K datasets from DeLaney et al. 
('2010') . The thick disk is apparent in the structure of the 
SNR as seen from the positive N axis (middle panel of 
Fig.^, and (to a lesser extent) the positive E axis (right 
panelof Fig. [t]). Among our observed LEs, LE2116 and 
LE3923 are sampling LoS away from this disk, and do 
not intersect any prominent ejecta structures. The LoS 
of LE2521, on the other hand, intersects the edge of the 
large complex of Fe K emission in the NW that is at the 
edge of the thick disk. 

6. DISCUSSION 

6.1. Systematic Effects 

To properly assess differences in the LE spectra, we 
have compared each LE to light-curve and window- 
function weighted spectra of well- studied SNe. When as- 
sessing the similarities of a LE to another SN, this process 
is straightforward. However, when comparing different 




Fig. 6. — Observed LE spectra near the He I A5876 (left), Ho; (middle), and Ca II NIR triplet (right) features, shown on a velocity scale. 
The red, blue, and orange curves correspond to LEs LE2116, LE2521, and LE3923, respectively (see Fig. [s] for details). The top LE2521 
spectrum (which is not shown in the right-hand panel) is from the MMT; the remaining spectra are from Keck. The solid and dashed black 
curves are the corresponding light-curve weighted integrated spectra of SNe 1993J and 2003bg, respectively. 
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Fig. 7. — Three-dimensional [Ar II] (red) and Fe K (blue) datasets from'DeLaney et al." ('2010'). The datasets are viewed from Earth's 
perspective (as they appear on the plane of the sky, left panel), from above (from a point along the positive N axis, middle panel), and 
sideways (from a point along the positive E axis, right panel). The datasets are represented in discrete voxels (volumetric pixels), roughly 
0.2^^ on a side (see DeLaney et al. 2010 for details). The expansion center is represented by a black cross in all three panels. The LE 
directions are shown for the fiducial distance of 3.4 kpc (solid lines), and for reasonable upper and lower limits, respectively 4.0 and 3.0 
kpc (dotted lines). 
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LEs by examining their similarities to other SNe, there 
are potential systematic effects. 

6.1.1. Window Functions and Comparison Objects 

When making integrated comparison spectra, the light 
curve and spectra of that particular SN are used. If that 
object is not a perfect analog of Cas A, then different 
window functions can introduce apparent differences be- 
tween the LEs. As an extreme example, consider two 
window functions: one which is a top hat only covering 
the premaximum portion of the light curve, and another 
top hat only covering the postmaximum portion of the 
light curve. If the comparison SN had lower line veloci- 
ties than Cas A before maximum brightness, but similar 
line velocities after maximum, the comparison spectra 
would indicate a difference for the two LEs. 

For the LEs presented here, the window functions are 
not significantly different. In Figure [3j we see that the 
major differences are the contribution of the light curve 
50-140 days after maximum brightness. One main source 
of uncertainty in the window function is caused by the 
uncertainty in the dust inclination. The dust inclination 
is derived from the apparent motion. In particular, for 
faint LEs, the uncertainty in the apparent motion can 
introduce u ncertainties in the inclination of ~ 5° (Rest 
et al.|[2QTQ ). We have not yet fully explored all ways to 
improve how we determine the dust inclination, and we 
hope that better modeling will further reduce its uncer- 
tainty in the future. Another source of uncertainty is 
caused by imperfections in the light-echo profile model- 
ing and thus the dust width caused by substructure in 
the dust. This can be tested with the SN 1987 A LEs 
dKunkel et al.||1987t |Crotts||1988| |Suntzeff et al.|[T988), 



since the spectrophotometric evolution of the SN itself 
has be en m onitored exten s ively. We show in [Rest et al . 
( |2010| and [Sinnott et al.| ( |201l| ) that the observed LE 
profile, apparent motion, and spectra can all be brought 
in excellent accordance with what is predicted using the 
SN 1987A spectrophotometric library. 

Figure [3] shows how the window functions change the 
integrated spectra, including a flat window function, cor- 
responding to the full light-curve weighted integrated 
spectrum. The most significant differences occur at the 
positions of nebular emission lines, particularly [O I] 
AA6300, 6363, [Ca II] AA7291, 7324, and the Ca II NIR 
triplet. 

The Ha line is most affected by the [O I] feature. For 
SN 1993 J, as [O I] becomes more prominent, it can make 
the velocity of the minimum of the Ha absorption appear 
lower, which may affect our interpretation of this feature 
for Cas A. However, it does not affect the blue edge of 
the emission feature. For SN 2003bg, the Ha absorption 
is at a higher velocity, coinciding almost exactly with the 
peak of the [O I] emission. As a result, the differences in 
window functions do not dramatically change the appar- 
ent velocity of the Ha absorption, but rather just affect 
its strength. 

For both SNe 1993J and 2003bg, the He I feature does 
not change much with different window functions. There 
is a noticeable change in the spectra between those with 
the measured window functions and the full light-curve 
weighted spectrum. The full light-curve weighted spec- 
trum has a slightly lower velocity for both SNe 1993 J 
and 2003bg. 



Although we cannot rule out that the window func- 
tion has created the differences seen in Cas A, that sce- 
nario requires particular tuning. One way this could hap- 
pen is if Cas A changes from a SN 1993J-like spectrum 
to a higher-velocity spectrum from about 80 days after 
maximum, when the window functions for LE2116 and 
LE3923 are 6 mag below peak, significantly less than the 
window function of LE2521. In addition, most of the flux 
at these late phases needs to be in Ha and He I instead 
of the nebular emission lines, which is physically unlikely 
for SNe lib. Then the integrated spectrum of SN 1993J 
could show a lower velocity than LE2521 while still hav- 
ing the same velocity for the other directions. Additional 
comparison objects and LE spectra with different window 
functions should address this potential systematic effect. 

6.2. Connection to Explosion Models 

The LE data alone suggest that Cas A was an asym- 
metric SN. As outlined in ^ this is consistent with ob- 
servations of the SNR, which indicate the existence of 



a rich set of multiple outflows (e .^., Hughes et al.||2000 



Willingale et al.||2002] [Fesen et al.||2006 a). The most 
prominent are bipolar Fe-rich outflows m the SE and NW 
that are blueshifted and redshifted, respectively. These 
outflows are approximately (but not exactly) aligned 
along a single axis through the kinematic center of the 



SNR jBurrows et al.|20Q5| [Wheeler et al.|2 008 DeLaney 
et al.lzOlol). The combination of LE and S 



y2010p . The combination of LE and SNK data paint 
a coherent picture of the SN, and further connect obser- 
vations to the explosion. 

Our detection of blueshifted Ha and He I emission from 
the direction of the NW outflow (LE2521) is the first di- 
rect connection between significantly higher kinetic en- 
ergy per unit mass in the SN explosion in one direction 
and outflows observed in the SNR. Additionally, the NW 
outflow that seems to be associated with LE2521 is in 
roughly the opposite direction of the apparent motion of 
a compact X-ray source, which is presumed to be the re- 
sulting neutron star (NS). From the position of the NS 
relative to the kinematic center of the SNR and the age 
of the SNR, Fesen et al.j ( |2006a| ) determined that the ob- 
ject is moving at 350 kms~^ in the plane of the sky with 
a position angle of 169° ± 8.4°. A simple explanation of 
both the high- velocity ejecta in the direction of LE2521 
and the NS kick direction and velocity is an inherently 
asymmetric explosion. 

Although this alignment could be coincidental, it may 
provide a key clue to understanding core-collapse SN ex- 
plosions. In particular, these observations indicate that 
the explosion mechanism may be directly connected to 
the observed SN — a nontrivial statement considering 
the size scales between the core and outer envelope of 
the star. 

One model to give the NS a kick is the eject a-driven 
mechanism, in which asymmetry in density before the 
collapse leads to an asymmetric explosion. This gives 
the protoneutron star a kick in the op posite direction 
than the ejecta, as in the rocket effect (Burrows & Hayes 
1996; F ryer 2004) . However, if neutrino-driven kicks help 
drive the supernova explosion, then the strongest ejecta 
motion can be in the direction of the NS, and opposite 



the direction of the neutrinos ( [Fryer fc Kusenko||2006 
"Socrates et al."2005>). 

Numerical models have shown that large asymmetric 
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instabilities c an occur near the core during collapse (e.g., 
Blondin et al. 2003; B urrows et al.|2007| [Marek &: Janka 



2009p . These instabilities can be imprinted on both the 



ejecta and NS by aligning high-velocity ejecta along 
the same axis as the NS kick. 3D models have only re- 
cently been explored, and the models have not been ex- 
tended to examine the effects of this mechanism on the 
composition, density, or velocity structure of the ejecta. 
Nonetheless, these mechanisms are an example of a pos- 
sible way to directly connect the NS kick with an asym- 
metric velocity in the ejecta. 

To determine if the alignment of the NS kick and 
higher-velocity ejecta can be physically explained, we 
perform a simple calculation using the ejecta-driven 
model. If we assume that in the direction of the NS's 
motion the ejecta velocity is the same as in the lower- 
velocity directions we have measured, and that the mo- 
mentum of the NS is equal to the excess momentum of 
the ejecta in the opposite direction, we can measure the 
mass of ejecta at higher velocity. The velocity difference 
is '^4000 kms~^ for the different LoS. We also assume 
that IM0 < Mns < 2M0. Then we have 



Ml 



HVej 



0.09 



/4000 kms"^ 



350kms"V V ^^hv 



Mns. (1) 



For our assumed NS mass, we find that the amount of 
material moving at the higher velocity would be 0.09 < 
Mnvej ^ 0.18 M0. This mass estimate assumes that 
the radial velocity of the NS is zero. Although the ra- 
dial velocity could be quite large, it is likely of order the 
transverse velocity, and will not affect our mass estimate 
by more than a factor of a few. Similarly, we have as- 
sumed that the velocity difference in the ejecta is the 
total velocity difference, but there could be a significant 
c omponent perpendicu lar to our line of sight. 



Willingale et al.|([2QQ 3) determined that the Cas A SNR 
had a total mass oi 2.2 with ~0.4 of fast-moving 
ejecta having an initial velocity of 15,000 kms~^, con- 
sistent with our measured velocity but slightly higher 
than our estimated mass of high- velocity material. Will-| 
ingale et al. (2003) claim that ~90% of the kinetic en- 



ergy of the SN was at the highest velocities, which seems 
inconsistent with SN explosion models. If the Cas A 
SN explosion was similar to these models, which expect 
most of the kinetic energy to be at lower velocities (e.g., 
~ 6000 kms~^), then the lower mass estimate for the 
high- velocity material seems reasonable. This implies a 
kinetic energy of ~2 x 10^^ erg located in high- velocity 
material. The estimate of the kinetic energy of Cas A is 
-(2-3) X 10^^ erg (iLaming fc Hwang|2003| for an ejected 
mass of ~2M0. 'rhese values are very similar to those 
derived for SN 2006aj, the SN Ic a ssociated with the X- 
ray flash GRB 06021 8 (e.g., |Pian et al.^ 2006) and which 



Mazzali et al. ( 2006[ ) suggested to have been the result 
ot a magnet ar event. The analogy with the properties 
of Cas A and the orientation of the motion of the NS 
suggests that the Cas A SN may have produced a mag- 
netar. This would justify the higher-than average kinetic 
energy of Cas A. 

A further constraint for explosion models is the differ- 
ent velocities of H, He, and Ca for LE2521 relative to 
SNe 1993J and 2003bg. SNe 1993J and 2003bg have 
similar velocities for the Ca II NIR triplet, and both 



are consistent with that of LE2521. The Ha veloc- 
ity of SN 2003bg is significantly larger than that of 
SN 1993J, and is consistent with that of LE2521. How- 
ever, SNe 1993 J and 2003bg have similar velocities for 
He I, with SN 2003bg having a slightly larger velocity. 
Conversely, LE2521 has a much larger He I velocity than 
that of either comparison object. 

There are several ways to explain the velocity struc- 
ture. First, since all spectra appear to have similar ve- 
locities for Ca II, it is likely that the feature is forming in 
a more symmetric region of the ejecta. The small differ- 
ences between the LEs and the relatively low velocities of 
the feature suggest that the Ca-emitting region is more 
central than that of H or He, and that the central region 
of the ejecta is more symmetric than the outer layers. 

There are several explanations for the He I velocity 
being higher than that of SN 2003bg, which had the same 
Ha velocity of LE2521. In Cas A the H layer may have 
been extremely thin (thinner than in either SN 1993J 
or SN 2003bg), causing the H and He velocities to be 
coincident. Alternatively, the Cas A He layer may be 
more mixed into the hydrogen layer than in SN 2003bg, 
causing He I to have a velocity similar to that of Ha. 
Finally, the ^^Ni distribution may have been different in 
the two objects, causing a different ionization structure 
in the outer ejecta. 

7. CONCLUSIONS 

We have obtained optical spectra of LEs from three 
different perspectives of the Cas A SN, effectively prob- 
ing different regions of the SN photosphere — the first 
time that this technique has been applied to a SN. The 
spectra are very similar to each other and are all similar 
to the prototypical SN lib 1993J. After accounting for 
the window function determined by the combination of 
dust inclination and slit orientation, we are able to pre- 
cisely compare Cas A to other SNe as well as compare 
the LEs to each other. From these comparisons, two 
of the three directions have spectra which are indistin- 
guishable from that of SN 1993 J; however, one direction 
has He I and Ha P-Cygni features that are significantly 
blueshifted (-4000 kms"^) relative to SN 1993J and the 
other two directions, indicating a higher ejecta velocity 
from Cas A in that one direction. This is direct and 
independent evidence of an asymmetric explosion. 

The spectrum for the discrepant LE has an Ha 
line profile consistent with that of the high-luminosity 
SN lib 2003bg, but its He I A5876 hue profile had an 
even higher velocity than that of SN 2003bg. This may 
indicate that Cas A had a very thin hydrogen layer, sig- 
nificant ejecta mixing, or different ionization structure 
in this direction. All LE spectra have Ca II NIR triplet 
line profiles consistent with each other as well as with 
those of SNe 1993J and 2003bg. This suggests that the 
emitting region of the Ca is distributed more spherically 
than that of the H or He-emitting regions. 

Even though there seems to be a "jet-like" structure 
in the NE corner and a co unterjet in the SW corner 
(Hwang et al. 2004 Fesen et al. 2006a[ ), recent optical 
and X-ray data from the Cas A feNK indicate that the 
dominant Cas A SN outflow is in the SE at a position an- 
gle of —115°, slightly tilted toward the observe r, and its 
count erpart approximately on the opposite side (Burr ows] 
et al.||2005| [Wheeler et al.^.2008, ,DeLaney et al., 20fUl r 
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Our detection of a blueshift looking into the counter 
outflow in the NW corner is the first direct, unambigu- 
ous, and independent confirmation of this outfiow. It is 
also in excellent agreement with the apparent motion of 
the compact object, which mo ves at a position angle of 
169° ±8.4° (Tananbaum 1 999| |Fesen et aL,2QQ6b) away 
from the center of the SNlT 

Finally, we note that the existing surveys for LE fea- 
tures in this portion of the Galactic plane are far from 
complete, and that additional LEs are very likely to be 
discovered, providing additional perspectives of the SN 
in three dimensions. The inventory of such features will 
further illuminate the degree of asymmetry of the SN, 
but will also serve the purpose of testing the degree of 
coherence of spectra from similar perspectives. Since the 
spectrum at a given dust-concentration location is the 
result of integration over an entire hemisphere of SN pho- 
tosphere, spectral differences are expected to vary slowly 
with changes in the perspective angle. 
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